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Abstract. We study the reformulation of integer linear programs by means of a mixed 
integer linear program with fewer integer variables. Such reformulations can be solved 
efficiently with mixed integer linear programming techniques. We exhibit examples that 
demonstrate how integer programs can be reformulated using far fewer integer variables. 
To this end, we introduce a generalization of total unimodularity called the affine TU- 
dimension of a matrix and study related theory and algorithms for determining the 
affine TU-dimension of a matrix. We also present bounds on the number of integer 
variables needed to represent certain integer hulls. 
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1 Introduction 

Reformulations of integer programs with linear constraints are common in the integer pro¬ 
gramming literature. The main motivation behind these reformulations is that linear programs 
can be rapidly solved in both theory and practice. Often integer programs are reformulated 
such that the linear programming relaxation is improved, hence creating better bounds for 
branch-and-bound based algorithms |22j . These reformulations include Lagrangian relaxation, 
Dantzig-Wolfe reformulation, and cutting planes. In a much stronger sense, many combinato¬ 
rial optimization problems can be formulated exactly using linear inequalities. These formu¬ 
lations typically involve totally unimodular (TU) matrices or totally dual integral systems. 
Since the feasible set of the linear relaxation exactly describes its integer hull in these settings, 
the optimization problem can be solved by solving the linear relaxation. See, for instance, m- 
The aim of this work is to instead reformulate integer linear programs by means of a mixed 
integer linear program with few integer variables. We focus on reformulating the feasible region 
of the integer linear program by mixed integer constraints. For a polyhedron P = {x G M" | 
Ax < &}, we wish to hnd a matrix W € such that 

conv(P n Z") = conv({a: £ P \ Wx £ Z^'}). (1) 

Since W is an integral matrix, the inclusion C in property (HD is always fulfilled. Furthermore, 
if W is the n x n identity matrix, the property (HD is trivially satished. The objective here is 
to find an integer matrix W with few rows k that can model the integer hull of the feasible 
region as in HD- 
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With such a reformulation at hand, the underlying optimization problem can then be 
solved using mixed integer linear programming techniques. From a theoretical point of view, 
Lenstra |14) presented an algorithm to solve mixed integer linear programs in polynomial time 
when the number of integer variables is fixed. Also in practice, mixed integer linear programs 
with few integer variables can lead to algorithms with improved running time m- 

As a centerpiece of this paper, we study the following generalization of total unimodularity 
that admits reformulations of P for all integral right hand sides b. 

Definition 1 We say a matrix A G admits a k-row affine TU decomposition if there 

is totally unimodular, where 

A G is the unique matrix satisfying A = A + UW. The minimum k such that A admits 

a k-row affine TU decomposition is called the affine TU-dimension of A. 

Observe that A has affine TU-dimension 0 if and only if A is TU. Also, the affine TU- 
dimension of any matrix is at most n since A = Omxn, U = A and W = In produces a valid 
affine TU decomposition. Affine TU decompositions with k rows can admit a model of the 
feasibility region of related polyhedra with only k integer variables. 

Theorem 2 Let A = A-\- UW G with W G {0, fee an affine TU decomposition, 

b G Z™ and £ G (ZU {—oo})", u G (ZU {oo})” with i < u. 

Then conv ({a; G R" | ^ < a; < u, Ax < b,Wx G Z^jj is an integral polyhedron. In particu¬ 
lar, P = {x G M" \ £ < X < u, Ax < 6} and W satisfy property ([T|). 

As an example of the convenience of Theorem [51 the parity polytope has a simple descrip¬ 
tion using one integer variable. 


exist matrices U G and W G {0, ±1}^^" such that 


Example 3 (Parity polytope) The n-dimensional (even) parity polytope is the convex hull 
of all 0-1 vectors in M" that have an even cardinality support. This polytope has exponentially 
many inequalities m, but can be described by an extended formulation with only 4n — 1 
inequalities mm- Alternatively, it can be described using one integrality constraint as P = 
conv ({a; G [0,1]" | | ^ ^}) • ^ projection of the polytope 


Q = conv 

= conv 


€ [0,1]" X 
G [0,1]” X ] 


Xi 2z = Q, z gZ 


i=l 


A 


<b,W 


G Z 


with A = 
position as 



2 

-2 


b = 


r 

0l 


■ 2 ■ 

n 

- 1 ^ 

L n 

1 

' O 

-f 

-2 


and W = The matrix A admits an affine TU decom- 

W. By Theorem\^ Q is an integral polyhedron, thus P is 


also an integral polyhedron. 


In Section [51 we expand upon the theory of affine TU decompositions: we investigate 
structure of affine TU decompositions and prove Theorem [2] and related results. 

In Section |3l we give various examples of how understanding the affine TU-dimension can 
create a mixed integer model with few integer variables. 

In Section |4l we focus on computational issues in connection to affine TU decompositions 
and study the complexity of determining the affine TU-dimension of a matrix. In particular, 
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we show that it is iVP-Hard to decide if the (affine) TU-dimension oi A G is less than 

n. When k and the number m of rows of A is fixed, we give a polynomial time algorithm to 
determine if the affine TU-dimension of A is equal to k. 

In Section [5l we study mixed integer reformulations for knapsack polytopes. We prove 
a general lower bound of ^ integrality constraints necessary to achieve property © for the 
linear relaxation of knapsack polytopes. We then give a nonconstructive proof that in every 
0-1 knapsack polytope we can replace the integrality on all the n variables by at most n — 2 
integrality constraints to achieve property ©. Apart from the added integrality constraints, 
we use only the original knapsack inequality and the 0-1 bounds on the variables, but no 
additional inequalities. In a final example we show the potential power of the addition of linear 
inequalities in the mixed integer model. We present a class of knapsack polytopes having an 
exponential sized polyhedral description. For these polytopes, by introducing linearly many 
additional linear inequality constraints one can replace the integrality constraints on all the 
variables by a single joint integrality constraint. 


Related Work 

The matrices with affine TU-dimension 1 have been called nearly totally unimodular matrices 
in [^. One of several applications of matrices with affine TU-dimension 1 is edge coloring in 
nearly bipartite graphs [6]. An undirected graph G is called nearly bipartite if it is not bipartite 
but one can obtain a bipartite graph by deleting one vertex of G [5] . The incidence matrix of 
a nearly bipartite graph has affine TU-dimension 1. By Theorem O integer hulls described by 
nearly totally unimodular matrices can be captured using just one integrality constraint, and 
hence integer linear optimization on these problems can be done in polynomial time using a 
mixed integer linear program. This, however, is not mentioned in [^. 

Mixed integer reformulations were studied by Martin m under the name variable redefini¬ 
tion. Martin showed how to reformulate problems that can be solved using dynamic programs. 
The main motivation here was to create tighter linear programming relaxations that improve 
bounds in a branch and bound algorithm. 

Many polynomial time algorithms in combinatorial optimization rely on the method of 
guessing the value of certain problem-related variables. Depending on the situation, guessing 
can for instance be done by polynomial enumeration of all possible combinations of values. 
Often one can interpret this approach as a polyhedral problem with an underlying mixed 
integer property like in ([T]). As an example, Hassin et al. [7] gave an efficient polynomial time 
approximation scheme for the constrained minimum spanning tree problem. Their algorithm 
relies on a partition of the edge set into (logarithmically many) buckets Ti. Then it guesses the 
number of chosen edges X^eeT- each bucket Ti, and constructs an approximate solution 
from this information. The guessing can be interpreted as solution of the linear relaxation of 
the polyhedral problem with additional integral constraints X^egT for all i. 

Another example of guessing integral values of certain variables is presented in Oriolo 
et al. m- They provide an exact efficient algorithm for a special case of a network design 
problem on rings which relies on guessing the right value for an integral variable b. Since their 
integral variable could have an exponential range, a complete enumeration is not efficient. 
They thus provide another argument, leading to an efficient procedure. The introduction of 
the integrality constraint & € Z gives an alternative conceptually simple way to obtain a 
polynomial time algorithm. 

We would like to mention the close connection of mixed integer reformulations to extended 
formulations. Consider the special case that the number of possible values d = Wx G in 
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© is polynomial in n. Then one can rewrite 

conv({a; £ P \ Wx G Z^'}) = conv | [J {a; G T* | Wx = d} 

\d&*‘ 

as a linear program of polynomial size in an extended space with a method of Balas [1]. In 
this special case, a formulation as in o can serve as a compact certificate for the existence of 
such an extended formulation. 

Lastly, we note that the set {a; G K" | Ax < 6, VFa; G is a projection from a n + fc 
dimensional space. When W is unimodular, as a convenient feature, we can find a change of 
variables that allows the integrality constraints Wx € to be modeled easily using integral 
variables without increasing the total number of variables. To see this, consider the Hermite 
normal form transformation of a unimodular matrix W G with full row rank k, that is 

W ■ L = [Ik Ofcx(ra-fc)] with L G Z"^" unimodular. Then by the bijective linear transformation 
X = Ly, we have 

maxjc^x I Ax < b, Wx G Z^} = max{(c^T)?/ | ALy < h, yi,... ,yk € Z}. 

Using this transformation prevents the need to write the problem in a higher dimension to 
model the integrality constraints. 

2 Properties of TU decompositions 

We will begin with a homogeneous version of affine TU decompositions that is easier to study. 

Definition 4 We say a matrix A G Z™^” admits a fc-row TU decomposition if there exist 
matriees U G Z™^^ and W G {0, ±1}^^" with A = UW such that W is totally unimodular. 
The minimum k such that A admits a k-row TU decomposition is called the TU-dimension of 

d. 

An obvious lower bound on the TU-dimension of a matrix is given by its rank, therefore even 
a TU matrix can have a large TU-dimension. The affine version of TU decompositions can 
rule out this artifact. 

Following Theorem [21 for a given matrix A, we are interested in hnding a fc-row affine TU 
decomposition A = A + UW with small k (the number of rows of lU). With this purpose 
in mind, U can be restricted to not contain an all-zero column. Otherwise, we could delete 
this r-th column of U together with the r-th row of W and still obtain a valid affine TU 
decomposition. 

Furthermore, we may assume without loss of generality that W G {0, has full rank 

k. This will allow us to simplify some proofs in what follows, but as the following remark 
shows it is not a real limitation. 

Remark 5 Let A G Z™^" and U G Z™^^ such that A = A + UW where 

W G {0,±1}^^" has rank k' < k. 

For any I C k}, let Wj. be the matrix eonsisting of the rows of W with indices from 

l, and U-i be the matrix eonsisting of the columns of U with indices from I. 

Now fix some I C |1,..., k}, with | J| = k' such that the TU matrix Wi. has rank k'. For 
any j € F := {1,... ,k} \ I, consider Wj., the j-th row of W. Then (Wj. = (Wi. Y" U G Z” 
for some U G . Since W is TU, there exists such a solution U with U G Z^ . 


A 

W 


is TU and 
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Then Wjc. = RWi. where R G fc )xfc is the matrix composed of rows (r^)^ forj G R. 

Since is TU, A = A + UW = A + U.iWi. + U.jcWic. = A+ (U.i + U.i.RfWj. is a 

W j. 

k'-row affine TU deeomposition of A. 


We now prove Theorem[2] A polyhedron P C R" is called integral if P = conv(P fl Z"). 

Proof (of Theorem\^. Let be the ri x n submatrix of the identity matrix with rows 
corresponding to the finite indices £i G Z, and let be the r2 x n submatrix of /„ with rows 
corresponding to the finite indices Ui G Z. 

Notice that any affine TU decomposition A = A-\-UW gives rise to an affine TU decomposition 


one has 


conv ({x G R" I ^ < X < u, Ax < b, Wx G Z^}) = conv ^|x G R" | Ax < b, Wx G Z^|^ . 



■ A ' 


■ A ' 


u 


b 

A := 

-R 

P 

— 

1 

to 

1 _ 

-h 

Opi xk 
pr2Xk_ 

W. Moreover, with b := 

1 

s 1 

1 _ 


Thus, it suffices to prove that for any affine TU decomposition A = A-\- UW, the set S = 
conv ({x G R" I Ax < 5, lUx G Z^}) is an integral polyhedron. 

Observe first that 5 is a polyhedron. To see this, let 

S = conv ({(x, d) G R" X Z^ I Ax < b, Wx — d = Oj,}) . 


Then S' is a polyhedron since it is the mixed integer hull of a rational polyhedron (see, e.g. (TU 
Section 16.7]). Since projections of polyhedra are also polyhedra, the projection S = proj 2 ,(S) 
of S onto the x variables is also a polyhedron. 

We now show that S is integral. Consider a decomposition of A as stated. For every fixed 
d G Z^, & — Pd is an integral vector, which implies that 


Pd ■= |x G R" I Ax <b — Ud, Wx = d| 


is an integral polyhedron by total unimodularity of 


'A' 

w 


[191 Theorem 19.1]. Then 


S = conv 



= conv 


[J conv(PdnZ") 


i.dGZ'' 


= conv 


U (Pd n z") 


I. 


= conv (S n Z"), 


and thus, S is an integral polyhedron. The second equality follows from the fact that Pd is an 
integral polyhedron and the last inequality follows since S' fl Z" = Ud^z^Pd H Z". This is clear 
since for any x G P fl Z", Wx G Z^ by integrality of W, and hence x G Pd for some d G Z^. 

Therefore, conv ({x G R" ] Ax < b, Wx G Z^'}) is an integral polyhedron. From the argu¬ 
ments outlined before, also conv ({x G R" \ (■ < x < u, Ax < 6, Wx G Z^}) is an integral poly¬ 
hedron. □ 

In a slightly restricted setting also a converse of Theorem [2] holds. 

For this, let us first link affine TU decompositions to affine decompositions containing 
general unimodular matrices. Recall that an integral (not necessarily square) r x s matrix 
is called unimodular if it has full row rank r and each nonsingular r x r submatrix has 
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determinant ±1. An integral matrix M G of rank r is unimodular if and only if the 

polyhedron {a; € | Ma; = r;, a; > 0} is integral for each v G Z™ (see, e.g. [TH Theorem 

19.2]). 


Now, let 


A 

W 


Im 

Ofcxm 


G ^(m+fc)x(ri+m) rank m + k. Notice that 


A 

W 


dm 

Ofcxm 


IS 


unimodular if and only if {(a:, y) G | Ax + y = b, Wx = d, a; > 0, y > 0} is integral for 

all b G Z™ and all d € Z^. The latter is true if and only if {a; G K." | Aa; < 6, ITa; = d, a; > 0} 
is integral for all b G Z™ and all d G Z^. 

For reference we sum up this insight as a remark. 


Remark 6 Let A = A + UW G be a decomposition of A with only integral matrices. 

\ A / I 

(i) The matrix ^ ^ ™ is unimodular if and only if {x G K" j Aa; < b, Wx = d, x > 0} 

is an integral polyhedron Jor all b G Z™ and all d G Z^. 

is unimodular then conv ({a; G M" | Aa; < b, Wx G Z^, x > O}) 


(ii) Moreover, if 


A 

W 


0fc> 


is an integral polyhedron for each b G Z™ . 

An anonymous reviewer pointed out to us that the equivalence in (i) above does not hold 
anymore if one drops the nonnegativity assumption on x. 

Theorem 7 Let A G and let W G have rank k such that the polyhedron 

{a; G K” I Aa; < b, Wx = d, x > 0} is integral for all b G Z™ and for all d G Z^. Then there 
exist matrices U G and W' G {0, ±1}^^" such that A = A + UW' is an affine TU 

decomposition. Moreover, for every b G ZM, 

conv ({x G K" I Ax < b, Wx G Z^, x > O}) = conv ({x G M” | Ax < b, W'x G Z^, x > O}) . 


Proof. Since W has rank k, we have n > k. By Remark [5] we know that 


IS a 


A Im 

kF O^cxm 

unimodular matrix. We claim that also W is unimodular. For this, consider any nonsingular 
k X k submatrix Wj of W, where J with | J| = fc is the set of columns chosen. Since U = 
Aj Ir 


Wj Ofcxm 
Let W = 


has rank m + k, we have det(Wj) = ±det(i3) G {±1}- Thus, W is unimodular. 

Wi W 2 ] where we assume without loss of generality that the last k columns of 
W are the full rank matrix W 2 . Then we have 

Ai A 2 I, 

Wi W 2 Ofcxm 


'A 



W 

O/cxm 





A 2 



Ofcxm 

W 2 



Ai - A2W2~^Wi 


Wi^Wi 


Omxk Im 
Ik ^kxm 


The matrix 


Im 

Ofcxm 


A 2 

LFz 


full row rank submatrix of the unimodular matrix 

T-l 


matrices, 


Ai - A2W2~"Wi 


is unimodular since it is (up to permutation of the columns) a 

. As a product of unimodular 
is unimodular as well. Moreover, it is well known 


'A 

W 


Im 

Ofcxm 


Omxfc Im 

W 2 ^Wi /fc Ofcxm 

that a matrix C G Z®^’’ is TU if and only if [C Is] is unimodular (see, e.g. [TU Section 19.1]) 

^mxk 


thus 


Ai - A2W2~"Wi 

W2~^Wi 


is TU. Trivially, the matrix 


Ai - A2W2~"Wi 


W2~^Wi 


is TU as well. 


In particular, A = [Ai —A 2 W 2 ^Wi 0mxfc]+A2[W2 ^Wi Ik] is an affine TU decomposition. 
Moreover, lUx G Z^ if and only if [W 2 "^kFi Ik\x gZ^. □ 
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3 Examples of TU decompositions 


Given a matrix, in order to find (affine) TU decompositions of it having few rows, one can 
exploit certain relations between its columns. We first show that any master knapsack polytope 

n 

t\iaster = conv({a; e {0,1}" I < b}) 

i=l 


can be modeled with only 0{^/n) many integrality constraints. Although there is a well- 
known polynomial time algorithm for integer linear optimization over PMaster using dynamic 
programming (see, e.g. [THl Section 18.5]), this characterization shows that the structure of 
the resulting problem is less complicated than that of general integer programs. 

Example 8 The affine TU-dimension of = [1,2,..., n] is 

It follows from Theorem\^ that PMaster can be modeled with 0{^/n) many integrality con¬ 
straints. 

Without loss of generality, suppose that n = for some I € Z_|_. This is without loss of 
generality since the affine TU-dimension of a"'" as a function of n is monotonically increasing 
and since 0{y/ri), 0{\y/n\), and 0(l^/n\) are eguivalent. 

We begin by exhibiting a matrix W to provide an upper bound. By defining 


and 


W = 


'h 

Oj 



h' 

Oj 


g{0^1}(2.-l)x^^ 


1 


T 

i J 


=[l 2 


2i ... e-i] 


one has a — l)-row TU decomposition = u^W. Thus has TU-dimension at most 
21 — 1. Since the first entry of the vector u is 1, we can set d to be the first row ofW, delete 
this row in W and the first entry in u, and get a {2£ — 2)-row affine TU decomposition. 

We next explain the lower bound. Consider any k-row TU decomposition = u^W with 
u € , and W G {0, ±1}^^" TU. Since the vector a has distinct entries, W needs to have 

distinct columns. As shown in Theorem 4.2], using slightly different terminology, a TU 
matrix with k rows has at most k'^ -\- k -\- 1 distinct columns (for a more modem approach, 
see also HE Section 21.3]). Since W has n distinct columns, k^ -\- k -\- 1 > n. In particular, 
k > y/n — 1. A simple argument shows that any k'-row affine TU decomposition of gives 
rise to a {k' -|- \)-row TU decomposition of . Hence, any affine TU decomposition for the 
master knapsack constraint vector = [1,2,..., n] needs to have at least ffn — 2 rows. 


Example 9 The affine TU-dimension of afi = [2, 2^,..., 2”] is n. 

Consider any TU decomposition of afi as = u^W where W has k rows and full row 
rank and u has no 0 entries. Since W is TU, it can be shown that there exist r^,... G 

{0, ±1}" that span the kernel ofW (see Lemma HE in Section^. By the decomposition above, 
= 0 for all i = 1,... ,n — k. But by the structure of a, observe that a^r 0 for all 
r G {0,±1}". Thus, we must have k = n. It follows that W is unimodular (and square) and 
that W~^ G Z"^". Thus we can write u = a. Since the greatest common divisor of the 

entries in is 2, u'^ has no ±1 entry. 




Now, consider any affine TU decomposition of as = of" + u^W and suppose the 

number of rows of W is less than n. Choosing and W = 

W provide a TU decomposition of a"’" with at most n rows in W. This is a contradiction with 
the above arguments since the first entry of is 1. Thus, the affine TU dimension of a"’" is 
at least n. Since W can be chosen as the identity matrix, the affine TU dimension of is 

exactly n. 


a 

W 


. But then and 


Next we give an affine TU decomposition for a class of matrices and show how this can be 
used to model the problem of finding r-flows in directed graphs. 

Example 10 (Block TU structure) 


(i) Let 


Ai 

be TU for i = 1,... ,r, where A'’ € with M = mi + ... + 

mr- Let W G for i = 1,... ,r and set 


A = 


A^ 


A^_ 

U’'A^ 


_U^A^ U^A2 . 

Then A has a (fci + ... + kf)-row affine TU decomposition with 

\A^ 


U = 


OMxfci 0Mxfe2 

C/1 t/2 


OMxfer 

Ur 


and W = 


A^ 


(ii) Consider the task of finding r flows in a directed graph G = (V,E) subject to certain 
joint capacity constraints as follows. Each flow /^: E —>■ Z (for i = 1,... ,r) is allowed 
to use the arcs Ei Cl E to satisfy a demand di to be transported from Si to ti. For each 
i we demand flow conservation with respect to fl at all vertices in V \ {si,ti}. Some 
of the arcs in E are allowed to be shared by several of the flow problems. For each arc 
e G E, define 1(e) = {i = 1,... ,r | e G Eflj to be the set of flow problems involving the 
arc e. We assume that for each arc e G E we are given a maximum capacity Ce € Z that 
acts as a joint bound /i(®) — ^e- 

This task is a variant of the multi-commodity flow problem. In contrast to our variant, 
in the multi-commodity flow problem one assumes that 1(e) = {1,..., r} for each e G E. 
As is standard for multi-commodity flow problems (see, e.g. m Chapter 19]), the above 
task can be modeled by using integral variables fi(e) G Z for eaeh e G E and each 
i = 1,... ,r. Linear constraints Af < b ensure feasibility of the flows and implement 
the capacity constraints. Notice that the set of capacity constraints decomposes into the 
ones for the arcs e G E with \I(e)\ = 1, and the ones for e with \I(e)\ > 2. 

One can rearrange the rows of A in order to get a matrix A as in part (i). For every 
i = 1,... ,r, the matrix A^ corresponds to the flow conservation constraints of the flow 
fi, together with the capacity constraints for e G E^ with |/(e)| = 1. Moreover, the U’’ are 
identity matrices, and the capacity constraints corresponding to e G E with |/(e)| > 2 
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form the matrix [A^ ... A^] where all the are subsets of the rows of identity matrices. 

lA^I 

Since is the constraint matrix for flow i, this matrix is TU. 


Applying the affine TU decomposition from part (i), we only need integral variables for 
the arcs that are allowed to be shared by multiple flows. These are X]eG-E-|/(e )|>2 — 

r ■ \{e G E \ |/(e)| > 2}| integral variables. 


Next we relate the class of almost totally unimodular matrices to the matrices with affine 
TU-dimension 1 (called nearly TU matrices by [5]). 

A square matrix A is called almost totally unimodular if A is not TU but every proper 
submatrix of A is TU. Almost totally unimodular matrices were introduced by Padberg [18] 
and are studied as a building block for fc-balanced matrices [U and references therein]. 

Example 11 Assume A = [a A] G is almost totally unimodular. Then one 

can write A = [0„ A] + a[l, 0,..., Oj. This is a 1-row affine TU decomposition, since every 
(n — 1) X (n — 1) submatrix of A is TU. Thus, every almost totally unimodular matrix has 
affine TU-dimension 1. 


This example shows that integer linear programs described with an almost totally uni¬ 
modular constraint matrix can be solved in polynomial time by solving a mixed integer linear 
program with only one integer variable. The authors of [23] study a more general class of 
matrices called almost unimodular matrices. Through studying the lattice width of polyhe- 
dra with constraint matrices that have small subdeterminants, they show that integer linear 
programs described by almost unimodular matrices can also be solved in polynomial time. 

We conclude this section with a discussion of the so-called integer decomposition property. 
A polyhedron P has the integer decomposition property, if for any positive integer k every 
integral vector in kP is the sum of k integral vectors from P. A matrix A G is TU 

(has affine TU-dimension 0) if and only if {a; G M" | Ax < b, x > 0„} has the integer 
decomposition property for each b G K"* [5]. It is straightforward from the proof in [5] that 
for a TU matrix A and b G lA, also P = {x G K." | Ax < 5} = conv({a: G IP \ Ax < 6}) has 
the integer decomposition property. Gijswijt [6] showed that for A G {0,±1}™^" with affine 
TU-dimension 1 and b G Z™, the polyhedron PA,b '■= conv({a: G Z" | Ax < 6}) has the integer 
decomposition property, too. 

The integer decomposition property does not hold in general for polytopes Pa,& where 
the matrix A has affine TU-dimension 2 and b G Z™. As a counterexample, consider the 
parity polytope P = conv ({(0, 0,0)^, (0,1,1)^, (1, 0,1)^, (1,1,0)^}). It does not have the 
integer decomposition property, since (1,1,1)^ G 2P fl cannot be written as the sum of 
two points in P n Z^ (see, e.g., 0)- The polytope P has a natural inequality description as 
P = {a: € I Ax < 6} with 


■ 1 

I 

1 ■ 


'2 

-I 

I 

I 

and b = 

0 

1 

-I 

I 

0 

1 

I 

-1 


0 


A has the 2-row affine TU decomposition 


■ 1 

1 

1 ■ 


■ 1 

0 

o ' 


■ I 

1 ■ 

-1 

1 

1 


-1 

0 

0 


I 

1 

1 

-1 

1 

— 

1 

0 

0 


-1 

1 

1 

1 

-I 


1 

0 

0 


I 

-1 
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Since P does not have the integer decomposition property, there exists no 1-row affine TU 
decomposition for A. Thus, A has affine TU-dimension 2. In particular, not every PA,b arising 
from a matrix A with affine TU-dimension 2 has the integer decomposition property. 

4 Determining the TU-dimension 

Seymour gave a decomposition procedure to recognize if a matrix is TU. A careful imple¬ 
mentation of this decomposition procedure results in an algorithm that decides in polynomial 
time if a matrix is TU, see Truemper m- Thus one can decide in polynomial time if a matrix 
has affine TU-dimension 0. Given any specific value k, we would like to decide if A has affine 
TU-dimension k. 

For B G let im(i?) := {By G R™ | y G R"} be the column space of B and 

ker(i?) = {x G R” | Bx = 0^} be the kernel of B. A TU decomposition of a matrix A G 
is strongly related to the existence of ±1 combinations of column vectors of A generating the 
vector Om, i.e. vectors r G ker(A) fl {0, ±1}" \ {0„}. 

We will need the following fact, a proof of which can for example be found in [Ill- 

Lemma 12 (Totally unimodular matrices have totally unimodular kernels) Let A G 

{0,±1}^” fc)xra with rankn—k. Then one can find a TU matrix W G {0, ±1}^^" 

such that ker(A) = im(IU^). 

We now show that it is AF-hard to decide if a matrix has affine TU-dimension n. 

Theorem 13 It is an NP-hard problem to decide if a matrix A G Z™^" has affine TU- 
dimension equal to n, even when we restrict to positive matrix entries and m = 1. 

Proof. It was shown in [21] that the equal-sum-subsets problem is AF-complete. Given b G Z", 
this problem is to decide if there exists a nonzero r G {0,±1}" such that b^r = 0. 

We reduce the equal-sum-subsets problem to deciding if the transpose of the vector 2n ■ b 
admits a (n — l)-row affine TU decomposition. In other words, for any vector b G Z", = 

2n ■ b^, we show that 6 is a yes-instance to the equal-sum-subset problem if and only if there 
exists a decomposition -I- u^W with u G W G {0, ±1}*^"“^^^", a G {0,±1}”, 

and W together with the row vector d^ is TU. 

If there exists an affine TU decomposition for as above with W TU having rank k < n—1, 
then by Lemma there exists a TU matrix R G {0, ±1}*^"“^^^” having rank n — k> 1 such 
that 

{xgR^\Wx = 0„_i} = {fU I t G R”-'=}. 

In particular R has only nonzero rows. Consider its first row r G {0, il}", which, as all other 
rows of W, satisfies Wr = 0„_i. Then a^r = a^r u^Wr = a^r, thus 

0 = o^r — aT’r = 2n ■ b'^r — r . 

Since d, r G {0, ±1}”, d^r G {—n ,..., n}. Furthermore 2n ■ b'^r G 2nZ, therefore aAr = 0 and 
b"^r = 0 for the nonzero r G {0, ±1}". 

Now assume there is a solution to the equal-sum-subsets problem, a nonzero r G {0,±1}" 
such that b"^r = 0. By Lemma 112] there exists a TU matrix W G {0, ±1}*-"“^^^" such that 

{x G R” I r'^x = 0} = {W'^u I u G R”-^}. 


11 


Thus b = W'^u has a solution u G and since W is TU we may choose u to be integral. 

Then = 2n ■ 6^ admits the (n — l)-row affine TU decomposition = 0^ + (2n • u)^lU, 
where W together with the row 0^ is TU. 

Therefore the problem to decide if a given matrix has affine TU-dimension n is fVP-hard. 

□ 


In view of the hardness result, it is an open question if there exists a polynomial time algo¬ 
rithm to decide if A has affine TU-dimension k when k is fixed. Below we provide polynomial 
time algorithms for other certain special cases of the task. We first present some auxiliary 
facts. A proof of the first one can be found in [m Section 19.3/19.4]. 


Lemma 14 For a nonsingular matrix E, 
is TU. 


B 

E 


D 

C 


is TU if and only if 


BE-^ 

-u-i 


D - BE-^C 

E-^C 


Lemma 15 Suppose A G admits an affine TU decomposition with W G {0, ±1}^^”. 

Then there exists a matrix W G {0,±1}^^" that contains Ik as a submatrix such that A also 
admits an affine TU decomposition with W'. 

Furthermore, there exists a unimodular matrix T G such that W = TW. The matri¬ 

ces T and lU' can he computed in polynomial time. 


Proof. Without loss of generality, assume that W has rank k. By reordering the columns of 
A and W with the same column permutation we may assume that Wi, the first k columns of 
W = [Wi W 2 ], have column rank k. 

Since Wi is unimodular, G is also unimodular. Then we have W := = 

[Jfe with = Wf G 

Consider an affine TU decomposition A = A-\- UW, where we define A = [Ai A 2 ] with 
Ai being the first k columns of A. In particular, the matrix 


Lemma [13] with E = Wi shows that 


matrix 


Omxfe A2 — AiW^ ^W2 

h 


AiWf^ 

_ -lUf 1 


A 2 - AiWf^W2 


Wf^W2 


W, W 2 ] 

is TU as well, and so is the 


Wf^W2 


. In particular, also W' is TU. Since we have 


A — [Ai A2] -I- U [lUi IU2] — [Omxfc ^2 — AiW;i ^W2] -I- [Ai UWi) \lk ^^ 2 ] , 

the latter is also an affine TU decomposition of A with the TU matrix W'. This completes 
the proof of existence for T := Wf^. The matrices T and W can be computed in polynomial 
time using basic linear algebra techniques. □ 


Theorem 16 Let n > k, Ai G , A 2 G and IU 2 G {0, TU. Then 


[Ai A 2 
only if 


admits a k-row affine TU decomposition with U G and W = [Ik W 2 ] if and 

A2 — A1IU2 


IU 2 


is TU. 


Moreover, given A and lU = [Ik W 2 ], we can find in polynomial time U and A such that 
A = A UW is an affine TU decomposition, if it exists. 


Proof. For the only if part, assume there is a matrix U such that [Ai 


A 

[Ik W 2 I 


is TU. Applying Lemma [Til on the TU matrix 


A 

[Ik W 2 I 


A 2 ] — A-\-U[Ik 
Ai -U A 2 
■ h 


W 2 ] and 
- UW 2 
tU2 
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with E = /fc implies that 


Ai -U A 2 - A 1 W 2 


is TU as well. In particular, 


A2 — Ai W2 
Wo 


-h W2 

is TU. 

For the if part, define A = [Omxfc A 2 — A 1 W 2 ] G and U = Ai & Then 


A 


Omxfc 

A2 — A1W2 

jik W2I 


. 

W2 


is TU by assumption, and [^1 A 2 ] = A + U[Ik W 2 ] by construction. 

Thus, in order to find an affine TU decomposition with given A and W = [Ik W 2 ] in poly¬ 
nomial time, we apply Truemper’s polynomial time algorithm ED to decide if 

is TU. If this is the case, we have an affine TU decomposition A = A + UW with A = 
[Omxfc A 2 — ^iIU 2 ] and U = Ai. Otherwise, there does not exist such a decomposition . □ 


A 2 — ^lW2 
IU2 


Observation 17 Assume that n is fixed and A G ig given. 

For an affine TU decomposition one has W G {0,±I}^'^" and in view of Remark[^ without 
loss of generality W has rank k < n. Moreover, by Lemma 1151 one can limit the search to 
matrices W that contain Ik as a submatrix. Enumerating these polynomially many matrices 
for all k < n, and applying for each of them the efficient algorithm from Theorem \16i one can 
find a matrix W having fewest number of rows such that A admits an affine TU decomposition 
together with W. 

Thus, one can compute the affine TU-dimension of A in polynomial time for fixed n. 


Next we show that for fixed k, the number of rows of W, and fixed m, the number of 
rows of A, we can decide efficiently if there exist U and W such that A admits an affine TU 
decomposition. 


Theorem 18 Suppose A G Z'”^". Then in polynomial time we can decide if A has affine 
TU-dimension k, provided that m and k are fixed. 


Proof. Let us first reduce the statement to the case that A has distinct columns. 

For this, without loss of generality assume the columns of A to be ordered as A = A 2 ], 

where Ai has distinct columns and A 2 has columns that are repeats of Ai. Let Ai have r 
columns. Consider any affine TU decomposition oi A = A-\- UW with a TU matrix [A W]. 
Let us divide the latter matrix into its first r columns, and the remaining ones 


'A' 


'ii 

A 2 

_W 


Wi 

W 2 


Then we can we can choose columns 


IT' 


as repeats of columns of 


Wi 


such that with 


Al = [Ai Ay and W' = [Wi Wy we get the (possibly different) affine TU decomposition 
A = A' -h UW'. In particular, A and Ai have the same affine TU-dimension. 

From now on assume that A G Z™^" has distinct columns and gives rise to an affine TU 

\ A 

decomposition A = A -\- UW. Then the TU matrix 


W 


with m -\- k rows also has distinct 


columns. A TU matrix with r rows and distinct columns has at most -|- r -|- 1 columns (for 

a reference, see Remark [5] above). Thus, n < (m -|- k)^ m k 1. 

Now with k and m, also n is constant. By Observation 1 171 we thus can compute the affine 
TU-dimension of A in polynomial time. □ 
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We close this section with a conjecture on the complexity of determining the affine TU- 
dimension of a matrix. 

Conjecture 19 Suppose A £ Then in polynomial time we can decide if A has affine 

TU-dimension k, provided that k is fixed. 

The above theorems and algorithms can be easily adapted to the consideration of TU 
decompositions instead of affine TU decompositions. In particular, the proof of Theorem [13] 
does also show the iVP-hardness of recognizing TU-dimension n. 

Observation 20 It is NP-hard to decide if a matrix A £ TU-dimension equal to 

n, even when we restrict to positive matrix entries and m = 1. 


5 Reformulations for knapsack polytopes 

In this section we investigate the size of reformulations specific to both A and b. In particular, 
we focus on bounds for the minimum number of integrality constraints needed to model a 
general 0-1 knapsack polytope. We then provide an example that demonstrates how adding 
even a single integrality constraint can vastly reduce the number of inequalities needed to 
describe the integer hull of a knapsack polytope. 

Lemma 21 (Lower bound) Let m £ Z_|_ and consider the knapsack polytope 


P = 


£ [ 0 , 1 ] 


2m 


2(xi -|- t/i) -|- 2‘^(x2 J/ 2 ) 2^{Xm + Vm) Si 2 


m+1 


-1 


in dimension n = 2m. Let W, W £ he integral matrices and let 


Q = conv ( P n 


[lU w'\ 


£ 


If k < m, then Q is not an integral polyhedron. 

Proof. Since k < m, dim(ker(IU)) > 1. Let r £ ker(IU) \ {0^} such that a^r > 0 where 


:= (2, 2^,..., 2™). Let x° £ {0,1}™ with x° = 


0 if Ti > 0 
1 if Ti < 0 


. Then there exists a A > 0 


with -I- Ar £ [0,1]™ for all 0 < A < A. Let be the complement of x°, that is, £ {0, l}'^ 

_ r2;0j 

and -I- = Im- Observe that — £ P fl {0, 1 }^™ since 


a'^x° + a^x^ = a'^lm = 2 ”"+^ - 2 < 2 "*+^ - 1 . 
Let 0 < A < minjA, where if a'^r = 0 we set -4^ to oo. 

Now, suppose that Q is an integral polytope. This implies that Qfl 


tegral polytope, and hence, its projection onto the first m variables, Q— := 


y = x^ 


X 

x^ 


> IS an m- 


e QK 


[0, Ij™ is also an integral polytope. By the choice of r and A, we have a^x'^ < ar{x^ + \r) and 

m 

a^{x'^ + Xr) + a^x^ = ^ 2* -|- Xa^r < 2™ — 1, 
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thus a;° + Ar S Q-^. Consider the optimizers G = argmaxja’^a; | x G Q^} and F = 
argmax{a^a; | x G {0,1}"*}. Since is integral, G = conv(F) 

By choice of a, the value a'^x is distinct for all x G {0,1}™ and by construction of Q^, 
F = {x°} and G = conv(F) = {x°}. But since a^x^ < a'^{x° + Xr), this implies that 
x'^ + Xr G G and hence + Ar = which is a contradiction with the fact that A > 0 and 
r Om- 

Thus, we conclude that Q was not integral. □ 

We next provide a positive result for modeling knapsack polytopes and show that every 
0-1 knapsack polytope can be modeled using at most n — 2 integrality constraints together 
with the upper and lower bounds on the variables. We first prove a theorem about separation 
of disjoint polytopes. 

Theorem 22 Let P,Q F R" he polytopes sueh that P fl Q = 0 and dim(conv(P U Q)) = n. 
Then there exist h G R", 01,02 G R with oi < 02 such that P G {x G R" | < Oi}, 

Q C {x S R" I h'^ X > 02}, and max{|{x € vert(P) | h'^x = Oi}|,|{x S vert((5) | h'^x = 
«2}|}> 

Proof. If either P or Q is empty, then the theorem is immediate. Assuming Q = 0, one may 
take a valid inequality h^x < oi for which {x G P \ = oi} is a facet of P, and then set 

02 > Oi. 

So suppose P,Q Let x G relint(P), y G rel int((5). Let P = P + t and Q = Q + t where 
t = X — 2y. Let x = x -I- t and y = y + t. Then x G relint(.P), y G relint((3) and x = 2y. We 
will prove the result for P and Q, which implies the result for P and Q. 

Since P fl Q = 0, by the (strict) hyperplane separation theorem, there exists a pair (h, a) 
such that h^X < d for all x G P and y > a for all y G Q. Define 

H = {{h,a) G R" X R I h^ x < a,h^y >o,o>d VxG vert(P),i/ G vert((3)}. 

Since (h, a) was a strict separator and since vert(P) and vert(Q) are finite, there exists an 
e > 0 sufficiently small such that the point (h, d -I- e) strictly satisfies all inequalities of H. 
Hence, {h,a -f e) G int(P), and therefore FI is full-dimensional. 

We claim that P is a bounded polyhedron. Indeed, let {f,P) G rec(P) where 

rec(P) = {(r,/3) G R” x R | r^x < 0,r^y > I3,j3>0 VxG vert(P),?/ G vert((3)} 

denotes the recession cone of H. By convexity, f^x < 0 and f^y > ^ > 0. Since x = 2y, 
this implies that f^x = f'^y = 0 and ,5 = 0. Since x G relint(P) and y G relint(Q), we have 
that f'^x = 0 for all x G P and f^y = 0 for all y G Q. Since dim(conv(P U Q)) = n, there 
exist G conv(P U Q) and a A > 0 such that Xf = z^ — z^. By definition of conv(P U Q), 
there exist x^, x^ G P, G Q, Mi, 1^2 G [0,1] such that z® = /x^x® -I- (1 — for * = 1,2. 
Combining this, we see that 


A|lr||2 = Ar^r = Ar'^(yiX^ -k + (1 - Mi)y^ + (1 - M2)y^) = 0, 

where the last equality comes from the fact that f^x = 0 for all x G P and f^y = 0 for all 
y G Q. Hence f = 0„ and rec(P) = {(0„,0)}. Therefore Ff is bounded (see, e.g., [HJ Section 
8.2 (5) (hi)]). 

Finally, since P is a bounded, full dimensional polytope, there exists a vertex (h, d) of P 
such that the inequality a > d is not tight. Since this vertex is defined by n-l-1 tight inequalities 
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from H that are not the inequality a> a, there exists a set of n + 1 points corresponding to 
the tight inequalities, each of them being a vertex of either P or Q. By pigeonhole principle, 
there must be at least such vertices in either P or Q. This completes the proof. □ 

Consider such a face F of P respectively Q spanned by at least [vertices. We want 
to find a TU matrix W € {0, such that Wx = d with the same d £ for all x € F. 

Below we prove that for all 0-1 knapsack polytopes this is possible for n > 4 with k = n — 2, 
where P = conv {a; £ {0,1}” | a'^x < b} and Q = conv({0,1}"\P). Thus, for all 0-1 knapsack 
polytopes, it suffices to introduce n — 2 integrality constraints to satisfy property (Ill- 

Theorem 23 Let n > 4, a £ Z", 6 £ Z and P = {x G [0,1]" | a^a: < 6}. There exists 
W £ sueh that (IT|) holds. 


Proof. Let S = P H {0, 1}" and = {0,1}" \ S. Without loss of generality we may assume 
S' ^ 0 and S'V 0. 

Since the polytopes conv(S) and conv(S''^) are disjoint and conv(S U S°) = [0,1]", by 
Theorem there exists h £ R", 01,02 with oi < 02 such that conv(S) C {a; £ M" | 
h'^X < Oi}, conv(S'^) C {x £ R." | > 02 }, and |{x £ S | /i^x = Oi}| > or 

i{x £ s^ I hPx = 02 }! > ■ 

Since n > 4, > 3. In case that |{x £ S | hPx = oi}| > let x^,x^,x^ £ {x £ 

S I hPX = 01 } be disjoint vertices of S. Otherwise, choose x^,x^,x^ £ {x £ S° | h^x = 02 } 
to be disjoint vertices of S^. Consider 


A = 


(x^ — x^) 
(x^ — x^) 


£ { 0 ,± 1 } 2 >^". 


Since x^, x^, x^ £ {0,1}”, the columns of A are from the set {[0, 0]^, [0, ±1]^, [±1, 0]^, ±[1,1]^}. 
Hence, A is TU. Furthermore, the affine hull aff(x^,x^,x^) satisfies 

aff(x^, x^, x^) = {x^ -I- A^y | y £ R^} = {x^ -I- x | x £ R", Wx = 0„_2}, 


where the rows of W span the kernel of the TU matrix A. By Lemma [T^] we may choose 
W £ to be TU. Since W it TU, the set <52 = {x £ [0,1]” | Wx = z} is an integral 

polytope for every 2 ; £ Z"“^. By construction of W, for every z £ Z"“^ there exists an a £ R 
such that (52 C {x £ R" I h'^x = a}. Thus either C conv(S') or C conv(S''^). Therefore 
Pr\Qz = conv(P n Z") n (52 and 


conv(P n Z”) 



conv({x £ P I Wx £ Z" ^}). 


□ 


It is an interesting open question to determine the minimum number of integrality con¬ 
straints needed to describe any knapsack polytope. We now know this number lies between ttL 
and n — 2. 

Allowing the addition of polynomially many linear inequalities to the mixed integer refor¬ 
mulation is natural and likely advantageous in many scenarios. In this vein, the parity polytope 
in Section [ 1 ] has been formulated using a single integrality condition together with the bounds 
on the variables and one additional inequality. We conclude this section with a family of 
knapsack examples in dimension n, where the polyhedral description consists of OirP) linear 
constraints. Here, /c is a parameter controlling the weight distribution of the instances. Each 
of these knapsack polytopes can be described by only Oin) linear constraints together with 
one integrality constraint. 
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Example 24 (Knapsack with big and small weights) Let nGZ, kGli, k>3 and 
b G M+. Consider weights ai € M+ where < cLi < ^ for all i G S = s} and 

j^b < Qi < b for all i G B = {s + 1,... ,n}. To simplify notation, assume that < Oi < 

02 < • ■ • < Os < ^ and s > k. We consider P = G [0,1]” | af"x < 6}. 

One can show that the facet description of the polytope conv(P fl Z") is given by 

X > 0 

jGB 

s 

Xj + Xj + (fc — 1) • Xj < k for aU l = l,..., s — k + 1 

j=i j^B: j&B 

aj ^b—ai 

X! + X! + (1^1 “ ■ X! < |i?| for allRQS,l<\R\<k-l, 

jeR j€B: j&B 

aj 

where i (i?) is the index of an item from R having smallest weight Oj. 

It follows easily that conv(P fl Z") is the convex hull of all x G [0,1]" that satisfy 

S 

Xj < k for all i = 1,..., s 

j=i j^B\ j^B 

aj '>b—ai 

Xj e {0,1}. 

j&B 

Notice that the latter mixed integer description in fact is an affine TU decomposition. 
However, the above polyhedral description for the considered knapsack polytope is specific for 
the right-hand-side b and the fixed k. 
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